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INVESTIGATION OF HALO SATELLITE ORBIT CONTROL 


It is desired to maintain the communication station in a "Halo" 
orbit near the translunar llbration point of the Earth-Moon system. 

Such a Halo orbit is to be “visible" at all times from Earth, and (ex- 
cept for perturbations due to solar gravity and mean Barth-Moon orbital 
eccentricity) is to be a possible periodic orbit for an unpowered ve- 
hicle. A one-parameter family of such periodic orbits, which are all 
unstable, has previously been computed analytically by a truncated 
Linstedt method [Ref. l], and is shown in Fig. 1. 

The accuracy of the truncated analytic description of a typical 
"nominaf* Halo orbit is limited, leading to an acceleration error averag- 
ing about 10 fi g’ s, this being therefore the "cost" of a very tight 
control to the nominal path. 

Since a looser control would be quite satisfactory, the station- 
keeping problem is posed in the following way: 
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6r is the deviation from the nominal path, e{ t) is the known accelera- 
tion error along the nominal path, and u is the control acceleration; 
find u(t) to minimize 
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where the positive scalar k is a measure of the tightness of the 
desired control. The optimal three-axis control n takes the form 

u = -C(t)x - ?(t) , 

3*6 _3xl 

where C(t), b(t) are periodic quantities, dependent on k, and 
computable with the aid of "spectral factorization" [2], The resulting 
rms |6rj and |u| are plotted in Fig. 2, for single-axis as well as 

3-axis control, together with an "error settling time" (= time for ini- 
tial } 6 r^ | to decrease by a factor 10) which, like the rms |6r|, increases 

as k decreases. The squares in Fig. 2 indicate that the average 
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control can be reduced to less than 10 g's with an error settling 
time of less than three months. 

More important, the resulting limiting motion in the 3-axis control 
case provides an improved nominal, correctable as in Ref. 1 for the 
non-coramensurable perturbations, which permit an average control of 
only 10 g' s with much tighter control, corresponding to settling 
times of the order of 1 day. 

Lastly, the station-keeping cost will be slightly increased by 
imperfect knowledge of the state x. If it is assumed that a continuous 
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control u is accompanied by a small random acceleration error of known 
statistics, and that observation errors have known statistics, the control 
law must make use of a Kalman estimator 

* = *<0* + [|( t )+u] + - ■(*)*] > 

z being the usual linearized measurement vector, and, if H(t) is 
periodic, K(t) approaches a computable periodic gain matrix as t 
increases. The increase in the rms ju| is then computable. 





In the case of range measurements from Earth, the information matrix 
T 

H (t)H(t), averaged over a day, is essentially constant, with informa- 
tion perpendicular to the Earth-Moon direction provided by simultaneous 
measurement from stations with different latitude and longitude. In the 
case of range-rate measurements, the Earth's rotation provides significant 
[3] but generally different NS and EW position information. The inclina- 
tion, however, of the Moon's orbit to the equator spoils the periodicity 
of the information relative to the Halo orbit, which has a period of 
about two weeks not exactly commensurable with the Moon's orbital motion, 

A roughly comparable, but strictly periodic, information pattern is 

obtained if the inclination of the Moon's orbit is taken as zero, Assum- 
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ing a random control execution error with rms value 10 g' s and a 
correlation time of one day, as well as typical range and range- rate 
measurement errors, the increase in the rms ju| was not significant, 
even when reliance was placed entirely on range-rate measurement. 

A comprehensive account of this work was written by M. Ratner as 0 
Pb.D. dissertation [4], 
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